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Weakly Interacting Massive Particles (WIMPs) are one of the leading candidates for Dark Matter. 
For understanding the properties of WIMPs and identifying them among new particles produced 
at colliders (hopefully in the near future), determinations of their mass and their couplings on 
nucleons from direct Dark Matter detection experiments are essential. Based on our method for 
determining the WIMP mass model-independently from experimental data, we present a way 
to also estimate the spin-independent (SI) WIMP-nucIeon coupling by using measured recoil 
energies directly. This method is independent of the as yet unknown velocity distribution of halo 
WIMPs. In spite of the uncertainty of the local WIMP density (of a factor of ^ 2), at least an upper 
limit on the SI WIMP-nucIeon coupling could be given, once two (or more) experiments with 
different target nuclei obtain positive signals. In a background-free environment, for a WIMP 
mass of 100 GeV its SI coupling on nucleons could in principle be estimated with a statistical 
error of only ^15% with just 50 events from each experiment. 
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1. Introduction 



There is strong evidence that more than 80% of all matter in the Universe is dark (i.e., interacts 
at most very weakly with electromagnetic radiation and ordinary matter). The dominant component 
of this cosmological Dark Matter must be due to some yet to be discovered, non-baryonic particles. 
Weakly Interacting Massive Particles (WIMPs) x with masses roughly between 10 GeV and a few 
TeV are one of the leading candidates for Dark Matter (for reviews, see Refs. 

Currently, the most promising method to detect many different WIMP candidates is the direct 
detection of the recoil energy deposited in a low-background laboratory detector by elastic scat- 
tering of ambient WIMPs on the target nuclei [^. The differential rate for elastic WIMP-nucleus 
scattering is given by ^ : 
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(1.1) 



Here R is the direct detection event rate, i.e., the number of events per unit time and unit mass 
of detector material, Q is the energy deposited in the detector, F{Q) is the elastic nuclear form 
factor, /i(v) is the one-dimensional velocity distribution function of the WIMPs impinging on 
the detector, v is the absolute value of the WIMP velocity in the laboratory frame. The constant 
coefficient .s^ is defined as 
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(1.2) 



where po is the WIMP density near the Earth and Oq is the total cross section ignoring the form 
factor suppression. The reduced mass /Mr n is defined by 



niy +mN 



(1.3) 



where nij^ is the WIMP mass and mN that of the target nucleus. Finally, Vmin = 0^^/Q is the minimal 
incoming velocity of incident WIMPs that can deposit the energy Q in the detector with 
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(1.4) 



and Vmax is relared to the escape velocity from our Galaxy at the position of the Solar system. 

2. Estimating the SI WIMP-nucleon coupling 

Based on our work on the reconstruction of the (moments of the) velocity distribution of halo 
WIMPs [U, the integral over the one-dimensional WIMP velocity distribution on the right-hand 
side of Eq.dO]), which is the minus-first moment of this distribution, can be estimated by [Q] 
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Here v{Q) = a^/Q, 2(mm.max) the minimal and maximal (cut-off) energies of the experimental 
data set, respectively, r{Q^^^ = {dR/dQ)Q=Q^^^, and 4(2min, 2raax) can be estimated by 
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(2.2) 



2 



Constraining the SI WIMP-nucleon coupling from direct DM detection data 



Chung-Lin Shan 



where the sum runs over all events in the data set that satisfy Qa G [Gminj 2max]-' On the other hand, 
using the assumption that the spin-independent (SI) WIMP scattering cross section is the same for 



both protons and neutrons, the "pointlike" cross section Gq of Eq.(|1.2|) can be written as 



(2.3) 



where /p is the effective XXPP four-point coupling, and A is the number of nucleons in the nucleus. 
Substituting Eqs.(2. 1) and ( p^ ) into Eq.(|l.l|), it can easily be found that 
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Note that the factor (§ appearing in the denominator is the exposure of the experiment, which is 
dimensionless in natural units. It relates the actual counting rate to the normalized rate of Eq.(|0]). 
The WIMP mass niyr^ on the right-hand side can be estimated by our method described in Ref. [Q] 
using data from two experiments. r(2min) and /q can be estimated from one of the two data sets 
used for determining or from a third experiment. Recall that, due to the degeneracy between 
the local WIMP density po and the WIMP-nucleus cross section Oq, one can not estimate each one 



of them without using some assumptions. Hence, by using Eq.([2.4[), one will need to accept an 
assumption of the local WIMP density po? 

In Figs. 1, we show the reconstructed spin-independent WIMP-nucleon coupling |/p|rec ^ 
function of the input WIMP mass WJ^ Following our work on determination of the WIMP mass 
[0], ^^Ge and ^^Si have been chosen as two targets for estimating ni;^. In order to avoid calculating 
correlations between and /q, a second independent data set with ^^Ge (left frame) or ^^Si (right 
frame) has been chosen for estimating /q. 

It can be seen in Figs. 1 that the reconstructed |/p|je(, values are underestimated for WIMP 
masses ^ 100 GeV. For Ge this deviation is larger than for Si. This systematic deviation is caused 
by the underestimate of /q. It is worse for heavier nuclei, where events with higher recoil energies 
can contribute more. This deviation of /q could be reduced by extending the maximal cut-off 
energy 2max to higher energy range since the kinematic maximal cut-off energy is larger for heavier 
WIMP masses. On the other hand, although we used the same event number for both experiments 
the statistical error on /q estimated with Si is larger than that with Ge. 

Nevertheless, it can be seen in Figs. 1 that, first, in spite of this systematic deviation, the true 
value of l/plrec always hes within the la error interval. Second, for a WIMP mass of 100 GeV, one 
could in principle already estimate the SI WIMP-nucleon couphng with a statistical uncertainty of 
only ~ 15% with just 50 events from each experiment. 



'Note that the generalized moments {V){v{Qjxiin),v(Qmax)) are independent of the local WIMP density, po, as 
well as of the WIMP-nucleus cross section, Oq. Moreover, one does not need to know /; (v) in order to determine its 



moments via a generalization of Eq.(2.1). Every term needed in this paper e.g., ^(Qmin) and /n(2mini 2max). can be 
estimated either from a functional form of the scattering spectrum or from experimental data (i.e., the measured recoil 
energies) directly. More details about estimating r{Q^\^), /n(2min> 2max), their statistical errors and the other formulae 
needed can be found in Refs. [^. 

^The most commonly used value for the local WIMP density is po ~ 0.3 GeV/cm-'. However, an uncertainty of a 
factor of ~ 2 has been usually adopted: po = 0.2 — 0.8 GeV/cm^. 
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Figure 1: The reconstructed spin-independent WIMP-nucleon coupling l/pl^c as a function of the input 
WIMP mass ni^ in- The open (red) squares indicate the input WIMP masses and the true values of the 
SI WIMP-nucleon couplings. The open (blue) circles and the (blue) crosses indicate the reconstructed SI 
WIMP-nucleon couplings and the la statistical errors. The theoretical predicted recoil spectrum for the 
shifted Maxwellian velocity distribution function ||l|], [|]] with Woods-Saxon elastic form factor [|l]] 
(vo — 220 km/s, vv = 231 km/s) have been used. The WIMP-nucleon cross section has been set to be 
10^*^ pb. ^^Ge and ^**Si have been chosen as two targets for estimating m^. A second (independent) data set 
with ^^Ge (left frame) or ^^Si (right frame) have been chosen as the third nucleus for estimating /q. Each 
experimental data set has 50 events under the maximal cut-off energy Qmax chosen as 100 keV. 
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Figure 2: The reconstructed spin-independent WIMP-nucleon coupling |/p j^g^, and the reconstructedWlMP 
mass nij^ jec on the cross section (coupling) v.s. WIMP mass plane. The open (red) squares indicate the input 
WIMP masses and the true values of the SI WIMP-nucleon couplings. The open (blue) circles indicate the 
reconstructed WIMP masses and the reconstructed SI couplings on nucleon. The horizontal and vertical 
soHd (blue) Hnes show the 1(7 statistical errors on m^j^c and l/pl^gc' respectively. Parameters as in Fig. 1. 
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Combining the estimate for |/p|^ with that for m^, in Figs. 2 we show the reconstructed SI 
coupUng l/plrec t^c reconstructed WIMP mass rec on the cross section (coupling) v.s. WIMP 
mass plane. We emphasize that by our methods described in Ref. [Q] and here, one can estimate 
m;^ and |/pp separately without any assumption for the WIMP velocity distribution. 

3. Conclusions 

In this paper we have extended our method for determining the WIMP mass [J] to estimate 
the spin-independent WIMP-nucleon coupUng from the elastic WIMP-nucleus scattering exper- 
iments. This method is independent of the velocity distribution of halo WIMPs as well as (prac- 
tically) of the as yet unknown WIMP mass. The only information needed is the measured recoil 
energies from at least two experiments with different target nuclei and the local Dark Matter density 
as the unique assumption. 

These information combined with the reconstructed WIMP mass will allow us not only to 
constrain the parameter space in different extensions of the Standard Model of particle physics, but 
also to identify WIMPs among new particles produced at colliders (hopefully in the near future). 
Once one is confident of this identification, one can use further collider measurements of the mass 
and couplings of WIMPs. Together with the reconstruction of the velocity distribution of halo 
WIMPs this will then yield a new determination of the local WIMP density. On the other hand, 
knowledge of the WIMP couplings will also permit prediction of the WIMP annihilation cross 
section. Together with information on the WIMP density, this will allow to predict the event rate 
in the indirect Dark Matter detection as well as to test our understanding of the early Universe. 
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